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Abstract 

In the context of three-flavor neutrino mixing, we present a thorough study of the phenomenolog- 
ical constraints apphcable to three observables sensitive to absolute neutrino masses: The effective 
neutrino mass in Tritium beta decay (m^); the effective Majorana neutrino mass in neutrinoless 
double beta decay {mpp); and the sum of neutrino masses in cosmology (S). We discuss the 
correlations among these variables which arise from the combination of all the available neutrino 
oscillation data, in both normal and inverse neutrino mass hierarchy. We set upper limits on 
by combining updated results from the Mainz and Troitsk experiments. We also consider the 
latest results on m^^ from the Heidelberg- Moscow experiment, both with and without the lower 
bound claimed by such experiment. We derive upper limits on S from an updated combination of 
data from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite and the 2 degrees Fields 
(2dF) Galaxy Redshifts Survey, with and without Lyman-a forest data from the Sloan Digital Sky 
Survey (SDSS), in models with a non-zero running of the spectral index of primordial inflationary 
perturbations. The results are discussed in terms of two-dimensional projections of the globally 
allowed region in the {mp, mpp, S) parameter space, which neatly show the relative impact of each 
data set. In particular, the (in)compatibility between S and m^^ constraints is highlighted for 
various combinations of data. We also briefly discuss how future neutrino data (both oscillatory 
and non-oscillatory) can further probe the currently allowed regions. 

PACS numbers: 14.60.Pq,23.40.-s,95.35.+d,98.80.-k 
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I. INTRODUCTION 



Atmospheric [11,1111,13, solar 0, i, 0, i, i, US , reactor QEIITI and accelerator 
neutrino experiments have convincingly established that neutrinos are massive and mixed 
TgL 17| . World neutrino data (with_a sing le controversial exception are consistent with 



a three-flavor mixing framework |19l . I20l |21L l22l l23l . 2^ . parameterized in terms of three 
neutrino masses (mi,m2, ms) and of three mixing angles (6'i2, 6^23, 6'i3) Jisl], plus a possible 
CP violating phase 6 l2| 



Neutrino oscillation experiments are sensitive to two independent squared mass difference, 



Sm"^ and Am^ (with ^ Am^), hereafter deflned as 



LOepe 





mi,m2,m3 



H ,±Am^ 



(1) 



where fixes the absolute neutrino mass scale, while the cases +Am'^ and — Am^ identify the 
so-called normal and inverted neutrino mass hierarchies, respectively. Neutrino oscillation 
data indicate that Sm"^ ~ 8 x 10"^ eV^ Q and Am^ ~ 2.4 x lO'^ eV^ 0, El. They also 
indicate that sin^^iz ~ 0.3 [ij, sin^^ss ^ 0.5 0, [13, and sin^ ^13 < few%B. However, 
they are currently unable to determine the mass hierarchy (±Am^) and the phase S |26l |. 
and are insensitive to the absolute mass parameter /i in Eq. (jT]). 

The absolute neutrino mass scale can be probed by non-oscillatory neutrino experiments. 
The most sensitive laboratory experiments to date have been focussed on tritium beta decay 
|2i[l,[3l| and on neutrinoless double beta decay ^[^[ll- 

Beta decay experiments probe 

the so-called effective electron neutrino mass mp [35^ . 



'^13'^12"^l + '^13'^12^2 



(2) 



where c?- = cos 6 ij and s?- = sin 6ij. Current experiments (Mainz [2^ and Troitsk |87l |) 



provide upper limits in the range m^j < few eV j25|, |3^ 

Neutrinoless double beta decay (Oi/2/3) experiments are instead sensitive to the so-called 
effective Majorana mass m^/? (if neutrinos are Majorana fermions) (sol. l4ol|. 



Ci3Ci2^i + Ci3S]^2'^2e ^ -|- S;^3m3e 



103 



(3) 



where 02 and 03 parameterize relative (and unknown) Majorana neutrino phases All 
0i/2/3 experiments place only upper bounds on (the most sensitive being in the eV range 
(3^). with the exception of the Heidelberg- Moscow experiment ||42i], which claims a positive 
(bu t highly debated js^]) 0z^2/9 signal, corresponding to m/3/3 in the sub-eV range at best fit 

El [3. 

Recently, astrophysical and cosmological observations have started to provide indirect 
upper limits on absolute neutrino masses (see, e.g., the reviews |20, |3l|, [i^ ) , competitive with 
those from laboratory experiments. In particular, the combined analysis of high-precision 
data from Cosmic Microwave Background (CMB) anisotropies and Large Scale Structures 



(LSS) has already reached a sensitivity of 0(eV) (see, e.g., j46|, |47|, |48| ) for the sum of the 
neutrino masses E, 

S = m\ + m2 + m3 . (4) 

We recall that the total neutrino energy density in our Universe, (where h is the Hubble 
constant normalized to = 100 km s~^ Mpc~^) is related to S by the well-known relation 
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Ojy/i^ = E/(93.2 eV) [25[, and plays an essential role in theories of structure formation. It 
can thus leave key signatures in LSS data (see, eg., jio^]) and, to a lesser extent, in CMB 
data (see, e.g^ |5fll])- Very recently, it has also been shown that accurate Lyman-a (Lya) 
forest data j5l|, taken at face value, can improve the current CMB+LSS constraints on E 
by a factor of ~ 3, with important consequences on absolute neutrino mass scenarios js^ . 
Given the potentially large systematics of Lya data jsj], we will conservatively perform and 
compare global cosmological fits both with and without Lya constraints. 

There is a vast literature about the theoretical and phenomenological links and correla- 
tions connecting [via Eqs. (0)^©] the 3z/ mass-mixing oscillation parameters ((5m^, Am^, 6ij) 
with the observables sensitive to absolute neutrino masses {mp^mpp.T,); see, e.g., the ex- 
tensive bibliography in the reviews 0, H H El. These interesting correlations have 



often been studied and plotted in terms of an auxiliary variable, such as the lightest neu 
trino mass mughtest; in particular, inany relevant papers have focussed on constraints in the 
(m/3/3, miightest) plane (see, e.g., [53|,|5^). However, since mughtest is not directly measured, 
it can actually be removed from data analyses, by charting the results just in terms of the 
three measurable quantities (m^,m^^,S).^ 

In this work we undertake a global phenomenological analysis of the constraints applicable 
to the observables {mp^mpji, E), by using up-to-date experimental data and state-of-the-art 
calculations for all the relevant laboratory and astrophysical quantities. To our knowledge, 
this is the first attempt to fitting world neutrino data, both oscillatory and non-oscillatory, in 
the "standard" astroparticle physics scenario emerged in the last few years (i.e., standard 3?/ 
mixing plus standard cosmology). The results are shown and discussed in terms of allowed 
regions in each of the three 2-dimensional projections of the parameter space (m^,m^/3, E), 
for both normal and inverted hierarchy. Such projections neatly show the interplay among 
different observables, as well as the impact of each data set. In particular, a tension is seen 
to emerge between current constraints on E and mpp. Our synthesis of world neutrino data 
appears also useful for estimating the sensitivity required by future experiments in order to 
solve, e.g., the mass hierarchy ambiguity, or to probe Majorana phases. 

The structure of the paper is as follows. In Sees. II, III, IV, and V we describe our 
treatment of the relevant input from oscillation, beta-decay, Oi^2/5, and cosmological data, 
respectively. In Sec. VI we discuss the bounds on (m^,m^^,E) coming from neutrino os- 
cillation data only, with emphasis on the correlations among couples of observables. In 
Sec. VII we combine the previous results with the upper bounds on (m^,m^^, E) coming 
from laboratory experiments and cosmological CMB-I-LSS data (without Lya forest data 
In Sec. VIII we implement a possible lower bound on based on the claim in 



and describe its effect on the global fit. Finally, in Sec. IX we discuss the impact of recent 
Lya data 51, S^l- We draw our conclusions in Sec. X. 



II. INPUT FROM u OSCILLATION DATA 



In this work, we perform an updated global analysis of oscillation neutrino data including: 
Atmospheric Super-Kamiokande (SK) data 0; solar data from the Chlorine (CI) Gallium 
(Ga) 0,0,3, SK 0, and Sudbury Neutrino Observatory (SNO) [l^ experiments; reactor 



^ Similar considerations apply to the auxiliary variable /i in Eq. being rn^^Yitcst = ~ Sm? /2 for normal 
hierarchy, and JTifig^tcst = ~ IS.m? for inverted hierarchy. 
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data from the KamLAND lll|, [131 and CHOOZ j28| experiments; and KEK-to-Kamioka 
(K2K) accelerator data 3 3- We refer the reader to previous works Jsi^ for details of the 



analysis. We include the latest experimental data presented at the Neutrino 2004 Conference 
j^, 0, 0]. In particular, we make a analysis of the available solar and KamLAND 
data, which are sensitive to the mass-mixing parameters ((5m^, sin^ 6'i2, sin^ 6'i3). We derive 
the likelihood of the Am^ parameter from a graphical reduction and combination of the 
corresponding likelihoods determined by the SK atmospheric ^ and K2K accelerator [l5l | 
experiments in the 2z/ limit, ^ marginalized with respect to 623."^ The probability distribution 
of Am^ is crucial to obtain upper bounds on ^^13 from the CHOOZ reactor experiment js^. 
The results of our analysis of neutrino oscillation data are embedded in a fi^ function, 

xlsci^m^, Am^ sin^ 6^12, sin^ 6*13) = 

Xsolar+KamL('^"^^ slu^ 9u, slu^ ^13) + 

Xsk+K2k(^"^^) + XcHooz(^"^^ sin^ ^13) , (5) 

which provides the likelihood for the four oscillation parameters ((5m^, Am^, sin^ 612, sin^ ^13) 
needed for the calculation of the {mp.misp, E) observables. 

Figure 1 shows the function ^x^^c = Xosc ~^i^Xosc from the global p oscillation data fit, 
marginalized with respect to each of its four arguments. The condition ^x^sc — provides 
n-a bounds on each mass-mixing parameter. At best fit we find 



V5rn2 ^ g j X 10"^ eV , (6) 
VAm2 ~ 4.9 X 10^2 eV , (7) 



and 



sin^ ^12 - 0.29 , (8) 

sin^ ^13 ~ 0.007 . (9) 

Notice, however, that the slight preference for nonzero sin^ Q\2 in Fig. 1 is not statistically 
significant. 



III. INPUT FROM TRITIUM /3-DECAY DATA 

Updated determinations of the effective electron neutrino mass squared have been 
recently presented js^ for the Mainz 3^| and Troitsk js^ tritium /3-decay experiments. The 



experimental values are consistent with zero within errors: 



= -1.2 ±3.0 eV^ (Mainz) , (10) 
= -2.3 ± 3.2 eV^ (Troitsk) , (11) 



^ At present, a full Zv analysis of SK+K2K data could only be jointly performed by the two experimental 
collaborations, since given relevant information is either unpublished or not reproducible outside the 
collaborations. In any case, iv effects are known to produce only negligible or minor changes in the Am^ 
likelihood from SK+K2K data (within current bounds on d\z)-, and thus they are not expected to alter 
the results of this work in any significant way. 

The angle ^23 is not relevant for the calculation of the (m/3, mpji, S) observables, see Eqs. ©-Q- 
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where errors are at la level, and systematic and statistical components have been added in 
quadrature. The function associated to the above /9-decay data can be simply defined 



as: 



= I '"^ I , (12) 

where x ^"^^ '^I3,x are the central values and errors on the right hand side of Eqs. fIlOp 
and for X = Mainz and X = Troitsk, respectively. By restricting the domain to the 
physical region (m^ > 0), the Ax^ function relevant for our analysis becomes 

^xl,ximp) = X%ximp) - Xax(O) , (13) 
providing the following upper bounds at 95% CL. (A^^ ^ 4), 

< 2.2 eV (Mainz) , (14) 

< 2.1 eV (Troitsk) , (15) 

in agreement with the upper bounds quoted by the experimental collaborations ;38j . 

It is matter of debate whether the two experiments have common systematics, e.g., re- 
sponsible for the negative central value of m| (see 0] for a recent discussion). In the 
absence of a clear indication in this sense, we assume that the two experimental results are 
independent, and simply combine them through 

,Mainz 

(m^) + Axl 

, Troitsk 

(m^) . (16) 

From the above definition, we get a combined upper limit at 95% C.L., 

m/3 < 1.8 eV (Mainz + Troitsk) , (17) 

which is less conservative than the 3 eV upper limit recommended in However, as 

we will see, upper limits on mp in the 2-3 eV range are, in any case, too weak to con- 
tribute significantly to the current global fit in the (m^, m^^, S) parameter space, so that 
"conservativeness" is not (yet) an issue in this context. 



IV. INPUT FROM GERMANIUM 0u2(3 DECAY DATA 

Neutrinoless double beta decay processes of the kind {Z, A) {Z + 2, A) + 2e~ have been 
searched in many experiments with different isotopes, yielding negative results (see jl^ Is^l 
for reviews). Recently, members of the Heidelberg- Moscow experiment [H^ have claimed 
the detection of a Oz/2/5 signal from the ™Ge isotope j43,|4^0|- '^^^ claimed signal would 
correspond to a decay half-life T in years (y) within the following 3a interval [^43. . ^4^] : 

logio(r/y) = 25.081°;^^ (3a) . (18) 



If the Oz/2/5 signal is entirely due to light Majorana neutrino masses, the half-life T is 
related to the rrifsp parameter in Eq. Q by the relation 



where rrip is the electron mass and Cmm is the nuclear matrix element for the considered 
isotope 1345. 

Unfortunately, theoretical uncertainties on Cmm are rather large (see e.g. (13] ), and their- 

m 



somewhat arbitrary — estimate is matter of debate (see j5J, |59|, |60| and refs. therein). Our 
approach to estimate a central value and an error for matrix element Cmm relevant to 
^^Ge experiments is the following. We take from the list of "acceptable" Cmm calculations 
(i.e., calculations with no subsequently recognized errors or biases) discussed in the 
"extremal" estimates: \og^o{CZ?,/j~^) = -13.85 ^1] and log^oiCmm = -12.88 [62] . 
Then we assume that their half-sum and their difference define, respectively, the central 
value and the 3a uncertainty (in log scale), namely: logiQ(C^^/y^^) = —13.36 and 3a = 
13.85 - 12.88 = 0.97. We thus take 

logio(C™™/y-^) = -13.36 ± 0.97 {3a) (20) 

which embeds a rather conservative error estimate — the 3a uncertainty being almost equal 
to an order of magnitude variation in Cmm, as also evaluated in 'H^. The above central 
value for Cmm is close to recent detailed theoretical calculations ^63 ]. 

From the previous equations, under the assumption of a positive Oi/2/3 signal in Ref. ji^ 
3], we derive that 

1 1 
logio("^/3/3/eV) = logio(me/eV) - -logio((^mm/y"^) - -logio(T/y) (21) 

= 5.71 + i(13.36 ± 0.97) - i (25.08toit) (22) 

= -0.23 ±0.53 (3a) , (23) 

where, in the last line, we have symmetrized the asymmetric range 25.08lo;24 as 25.23 ±0.39 
(according to the recipe suggested in jg^l) before adding the errors in quadrature (see [H^l for 
the statistical rationale of this procedure). Equation (PHj) allows us to define the Ax|^(m/3^) 
function associated to the claim in Ref. |43l. l44| . and provides the 3a range 0.17 < m^^ < 2.0 
(in eV). ^ 

The claim in 3] has been subject to strong criticism, especially after the first pub- 



lication j42| (see j3J] and refs. therein). Therefore, we will also consider the possibility that 
T = 00 is allowed (i.e., that there is no Oz/2/3 signal), in which case the experimental lower 
bound on disappears. Since one cannot really reprocess the data of 0, 3| under the 



assumption of no signal, nor combine such data with the negative results of other experi- 
ments with comparable sensitivity to ji^, the definition of a Ax|^(m/3/3) function in 
the absence of a Oz/2/5 signal is somewhat ambiguous and, in part, arbitrary. Our approach 
simply consists in stretching the lower error in Eq. (PHj) to infinity, so that the lower bound 
on 171/313 disappears, while the upper bound remains in the 0(eV) ballpark indicated by the 
most sensitive experiments to date [3^. It is worth noticing that none of the main results 
of our work depends on the precise definition of the upper bound on m/313. 

In conclusion, we adopt the following two possible Oz/2/5 inputs (and associated Ax^^ 
functions) for our global analysis: 

logio(m/3^/eV) = -0.23 ±0.18 (Oz/2/3 signal assumed) , (24) 



Our estimated 3a ran ge m ge € [0.17,2.0] eV overlaps — but does not coincide — with the range mpp S 
[0.1,0.9] eV quoted in |43l l44l | . since we use different estimates for the central value and uncertainties of 
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logio(m/3/3/eV) = -0.23^^^^ (Oz/2^ signal not assumed) , (25) 

where errors are at la level. Concerning the unknown Majorana phases 02 and 03 in Eq. (0), 
we simply assume that they are independent and uniformly distributed in the range [0,7r]. 

V. INPUT FROM COSMOLOGICAL DATA 

The neutrino contribution to the overall energy density of the universe can play a relevant 
role in large scale structure formation and leave key signatures in several cosmological data 
sets. More specifically, neutrinos suppress the growth of fluctuations on scales below the 
horizon when they become non relativistic. A massive neutrinos of a fraction of eV would 
therefore produce a significant suppression in the clustering on small cosmological scales 
(namely, for comoving wavenumber k ~ 0.05 h Mpc~^). 

To constrain S from cosmological data, we perform a likelihood analysis comparing the 
recent observations with a set of models with cosmological parameters sampled as follows: 
cold dark matter (cdm) density ^cdmh"^ ^ [0.05, 0.20] in steps of 0.01; baryon density flbh'^ e 
[0.015, 0.030] (motivated by Big Bang Nucleosynthesis) in steps of 0.001; a cosmological 
constant G [0.50, 0.96] in steps of 0.02; and neutrino density Q^h"^ e [0.001, 0.020] in 
steps of 0.002. We restrict our analysis to flat A-CDM models, Qtot = 1, and we add a 
conservative external prior on the age of the universe, to > 10 Gyrs. The use of a method 
based on a database of models instead of Markov Chains (see e.g. ji^) has the advantage 
of being more reliable in the definition of likelihood confidence contours at more than 3a 
level. No alternative model to a cosmological constant is considered: the dark energy is 
described by an unclustered fluid with constant equation of state {w = —1). Variations in 
the equation of state affect mainly curvature parameters (see, e.g., and therefore do 
not alter our results on S. 

The value of the Hubble constant in our database is not an independent parameter, 
since it is determined through the flatness condition. We adopt the conservative top-hat 
bound 0.50 < h < 0.90 and we also consider the la constraint on the Hubble parameter, 
h = 0.71 ±0.07, obtained from Hubble Space Telescope (HST) measurements ji^. We allow 
for a reionization of the intergalactic medium by varying the CMB photon optical depth Tc 
in the range Tc G [0.05, 0.30] in steps of 0.02. 

We restrict the analysis to adiabatic inflationary models with a negligible contribution 
of gravity waves. We do not consider the possibihty of isocurvature perturbations [g^] or 
topological defects (see, e.g., jH^). We let vary the spectral index n of scalar primordial 
fluctuations in the range n G [0.85,1.3] and its running dn/d\nk G [—0.40,0.2] assuming 
pivot scales at fco = 0.05 Mpc~^ and ko = 0.002 Mpc~^. We rescale the fluctuation amplitude 
by a prefactor Cuo, in units of the value C^^^^ measured by the Wilkinsin Microwave 
Anisotropy Probe (WMAP) satellite. Finally, concerning the neutrino parameters, we fix 
the number of neutrino species to N^, = 3, all with the same mass (the effect of mass 
differences compatible with neutrino oscillation being negligible in the current cosmological 
data jfi^). An higher number of neutrino species can weakly affect both CMB and LSS data 
(see, e.g., (t^) but is highly constrained by standard big bang nucleosynthesis and is not 
considered in this work, where we focus on 3z/ mixing. 

The cosmological data we considered comes from observation of CMB anisotropics and 
polarization, galaxy redshift surveys and luminosity distances of type la supernovae. For the 
CMB data we use the recent temperature and cross polarization results from the WMAP 
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satellite '46| using the method explained in [71] and the publicly available code on the 
LAMBDA web site j?^ . Given a theoretical temperature anisotropy and polarization angular 
power spectrum in our database, we can therefore associate a Xwmap to the corresponding 
theoretical model. 

We further include the latest results from the BOOMERanG-98 j?^, Degree Angular 
Scale Interferometer (DASI) IH, MAXIMA-1 Cosmic Background Imager (CBI) [t^, 
and Very Small Array Extended (VSAE) 0| experiments by using the publicly available 
correlation matrices, window functions and beam and absolute calibration errors. The CMB 
data analysis methods for the non-WMAP data have been already extensively described in 
various papers (see, e.g., [t^) and will not be reported here. 

In addition to the CMB data we also consider the real-space power spectrum of galaxies 
from either the 2 degrees Fields (2dF) Galaxy Redshifts Survey or the Sloan Di gita l Sky 
Survey (SDSS), using the data and window functions of the analysis of [t^ and ]47l|. We 
restrict the analysis to a range of scales over which the fluctuations are assumed to be in 
the linear regime {k < 0.2h~^ Mpc). When combining with the CMB data, we marginalize 
over a bias b for each data set considered as an additional free parameter. 

We also include information from the Lya Forest in the SDSS, using the results of the 
analysis of jH^ and jsj], which probe the amplitude of linear fluctuations at very small scales. 
For this data set, small-scale power spectra are computed at high redshifts and compared 
with the values presented in pi\. As in JH^I, we do not consider running. 

We finally incorporate constraints obtained from the SN-Ia luminosity measurements of 
IsO ] using the so-called GOLD data set. Luminosity distances at SN-Ia redshifts are com- 
puted for each model in our database and compared with the observed apparent bolometric 
SN-Ia luminosities. 

In Fig. 2 we plot the likelihood distribution for S from our joint analysis of CMB + SN- 
Ia + HST + LSS data, transformed into an equivalent Ax| function, which allows to derive 
bounds on S at any fixed confidence level. We take LSS data either from the SDSS or the 
2dF survey (dashed and solid curves, respectively).^ 

As we can see, these curves do not show evidence for a neutrino mass (the best fit being 
at S ~ 0) and provide the 2a bound S ^ 1 .4 eV. Such bound is in good agreement with 
previous results in similar analyses |46l . 47, 81, 8^, 13] ■ However our results can be considered 



complementary a nd, in some cases, independent from those analyses. For example, the recent 
SN-Ia data from [80| were not included in those analyses, so that a model with fi^ = 1 and 
break in the spectrum as in ^84] could produce a good fit to the CMB and LSS data, strongly 
increasing the upper limits on S. Models of this kind are ruled out by the inclusion of the 
SN-Ia data. 

We found that the CMB+LSS results in Fig. 2 are stable under the assumption of running, 
i.e. there is a weak correlation between the running and the neutrino masses in the range 
of scales probed by these data sets. For definiteness, we will use the combination including 
2dF data (solid line in Fig. 2) in the global analyses of Sees. VII and VIII. 

Also plotted in Fig. 2 is the Ax| function from a joint analysis of CMB + SN- 
Ia + HST + 2dF -|- Lya. No running is assumed in this analysis, and we find a 2a bound 
S < 0.47 eV, in very good agreement (despite the more approximate method we used) with 
the analysis already presented in [s^ . 



^ For the sake of brevity, the subdominant block of data (SN-Ia + HST) is not exphcitly indicated in figure 
labels. 
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As shown in Fig. 2 and already discussed in j52], the inclusion of the Lya data from the 
SDSS set greatly improves the constraints on S. However, we remark that the constraints on 
the linear density fluctuations obtained from this dataset are derived from measurement of 
the Lya flux power spectrum Ppik) after a long inversion process, which involves numerical 
simulations and marginalization over the several parameters of the Lya model. Since the 
effects of possible systematics need still to be explored further both observationally and 
theoretically, we here take a conservative approach, and discuss the implications of this 
results for neutrino physics in a separate section (Sec. IX). 



VI. CONSTRAINTS ON (m/3,m^/3,S) FROM OSCILLATION DATA ONLY 

In this Section we present and discuss the regions allowed by neutrino oscillation data 
in the parameter space (m/3,m^^,S). Figure 3 shows the projections of the 2a regions 
(AxQgj,=4), onto each of the three coordinate planes (m/3/3,S), {mpp.mp), and (m/3,S), for 
both normal hierarchy (thick solid curves) and inverted hierarchy (thin solid curves) . In each 
of the three coordinate planes, the observables appear to be strongly correlated. Such cor- 
relations can be largely understood in the approximation sin^ ^13 ~ 0, and by distinguishing 
the following three main cases for Eq. ([Q): a) ^ Am?] b) /i^ ^ Am^; and c) /i^ < Am?. 
Earlier discussions of correlations in the {mpj^, S) and {mpp^mj^) planes have been presented 



m 



isl 8^ and in js^, respectively. 



A. Degenerate spectrum (DS) 

For /i^ ^ Am^ in Eq. (^, the neutrino masses m^ {i = 1, 2, 3) form a degenerate spectrum 
(DS) and, in both hierarchies, for sfg ^ one has 

m ~ (/i,/i,/i) , (26) 

~ /i , (27) 

m^p ~ fi\cl^ + sl^e"^'\ , (28) 
S ^ 3/i . (29) 

By eliminating the auxiliary mass scale parameter /i, one gets the following linear correlations 
among observables: 

- f , (30) 

mf3f3 - f / ' (31) 
mpp ^ mpf , (32) 

where the factor 

/ = |c?2 + 4e^*^| , (33) 
which can take any value in the range c^2 — — / — 1) embeds our ignorance about the 
Majorana phase 02- 

Equation (j3Up explains the tight correlation in the (m^, S) plane of Fig. 3, where, for DS 
masses, the allowed region reduces to a "line." The correlation is instead relaxed by the 
variable factor / in the (m^/3, S) and (m/3/3, m^) planes, where the allowed regions appear as 
"strips" in the DS limit. 
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B. Partially degenerate (PD) spectrum 



For ^ Am^ in Eq. the neutrino mass spectrum can be considered as partially 
degenerate (PD), in the sense that the largest (smallest) mass splitting can(not) be resolved. 
For vanishing 5^3 and Sm^ one has 



(/i,/i, ^fi^±Am^) , (34) 
~ /i , (35) 

mpp ~ /^l4 + 4e**'l , (36) 



E ~ 2/i + ^/i2 ± Am2 , (37) 

where the upper (lower) sign refers to normal (inverted) hierarchy. By eliminating /x, the 
following correlations are obtained: 



S ~ 2m[3 + ± Am2 , (38) 
mpf^ c:^ nifsf , (39) 



2mf3,3 + Jmjp± fAm^ 
E ~ , (40) 

where / is defined as in Eq. fjHH|) . 

According to Eq. ()H8|1 . the regions allowed in the [mp, S) plane of Fig. 3 for normal and 
inverted hierarchies — which overlap in the degenerate spectrum case — branch out when the 
spectrum becomes partially degenerate and sensitive to ±Am^. In particular, the curve for 
inverted hierarchy bends upwards and ends at (m/j^min, ^mm) — (vAm^, 2VAm2), while the 
curve for normal hierarchy bends downwards, and eventually enters the regime of hierarchical 
spectrum discussed in the next section. The two hierarchies are instead not distinguishable 
in the {mf^p, mp) plane of Fig. 3, where Eq. ()39|1 provides, in the PD case, the same hierarchy- 
independent linear correlation as in the DS case [Eq. (j32p ]. 

Finally, Eq. (j4Uj) explains the correlation in the plane (m^^,S); in particular, by taking 
extremal values of / (1 and c^2 — s^g) at fixed mpp, one gets upper and lower bounds on E 
in the PD case, which are different in the normal hierarchy {+Am?) and inverted hierarchy 
[—Aw?) cases js^ H^; analogously, for fixed S one gets upper and lower bounds on m^p. 
In conclusion, it appears that for PD spectra the cases of normal and inverted hierarchy 
might be discriminated (the better the lower the neutrino masses) in the planes (m^, S) and 
(m^^, E), but not in the plane (m/3^,m^). 



C. Hierarchical spectrum (HS) 

Finally, only for normal mass hierarchy one can also consider cases with /x^ < Am?, 
where the smallest mass splitting cannot be neglected, and a truly "hierarchical spectrum" 
(HS) with three resolved masses (mi < m2 <^ ma) is obtained. For 5^3 ~ one gets in the 
HS case: 
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mp ~ ^^1^ - 5m2(cf2 - s?2)/2 , (42) 




Eliminating the auxiliary mass scale jj, from the above equations does not lead to particular 
transparent formulae, and we shall limit ourselves to a few comments. For /x^ = 5m^/2 

one gets from the above equations the minima m/3_min — \/sl^^m? and Smin — VSm^ + 

Am? + S'm?/2, but not the minimum of mpp. In fact, mpp^^i^ ~ is reached at e*'^^ = — 1 
and for values of /z^ slightly higher than 8rr? 12 in Eq. ()4H|1 . which also lead to values of mp 
and S above their minima. Therefore, while the allowed region for normal hierarchy has an 
endpoint at (m/j^min, Smin) in the (m^, S) plane, it continues indefinitely (logm^g/j — > — oo) in 
the other two planes. We also mention that, in the HS case, corrections induced by nonzero 
values of (not included in the above equation, but included in the fit of Fig. 3) can be 
nonnegligible, and contribute to the spread of the allowed regions for small values of the 
three parameters (m^,m^/3,E). 

We conclude this section with a few remarks on prospective improvements. Besides 
current neutrino oscillation experiments, a number of new experiments have been planned for 
the next decade, aiming at a better determination of the mass-mixing neutrino parameters 
[ssj l. and in particular of 6*13 [s^, which is currently unknown. A nonzero value of 6*13 is 
crucial to discriminate normal and inverted hierarchies (±Am^) through matter effects in 
future long-baseline experiments jo^]. If the neutrino oscillation parameters were all perfectly 
known, the allowed regions in the upper panel of Fig. 3 would shrink to two lines (for the two 
possible hierarchies). In the other two panels the allowed regions would be only moderately 
reduced, given the irreducible spread in m^p (for any fixed value of mp or S) induced by 
the unknown Majorana phases, and in particular by the phase 02 through the factor / in 
Eq. 

VII. ADDING UPPER BOUNDS TO (m/3,m/3^,S) FROM LABORATORY AND 
COSMOLOGY 

The results of the previous sections have been obtained by projecting Ax^ = 4 regions 
from neutrino oscillation data only (A^^ = Ax^^^). Here we consider also the regions ob- 
tained by projecting the various Ax^ functions associated with the mp data input (Sec. Ill), 
with the input in the conservative case [Sec. IV, Eq. (j^ ] and with the cosmological 
CMB+LSS input (Sec. V), both separately and in combination with Ax^^^.^ 

Figure 4 shows such projections in the three coordinate planes. Separate laboratory and 
cosmological upper bounds at the 2(j level are shown as dashed lines, while the regions 
allowed by the combination of laboratory, cosmological, and oscillation data are shown as 
thick solid curves for normal hierarchy and as thin solid curves for inverted hierarchy. It 
can be seen that the upper bounds on the (m^,m^^,E) observables are dominated by the 



^ We will add Lya data to the global fit in Sec. IX. 
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cosmological upper bound on S. This bound, via the (m/g,E) and S) correlations 

induced by oscillation data, provides upper limits also on and m^, which happen to be 
stronger than the current laboratory limits by a factor ~ 4. Since significant improvements 
on laboratorvlimits for m^^ and m^j will require new experiments and several years of 
data taking |38L l9l|. cosmological determinations of S, although indirect, will continue to 
provide, in the next future, the most sensitive upper limits (and hopefully a signal) for 
absolute neutrino mass observables. 

Figure 4 is useful to estimate the prospective impact of possible future measurements 
of the kind M ± (Tm, where M is any one of the three observables {mp^nijjp^Tj). The 
intersections of the currently allowed regions with the prospective M ± gm band(s) provide 
immediately an estimate of "future" allowed regions in the presence of one or more positive 
signals. In particular, we invite the reader to evaluate by herself or himself the relative 
accuracy needed to discriminate — if possible at all — the two mass spectrum hierarchies in 
the three planes of Fig. 4, as well as the m^js accuracy needed to reduce the vertical spread 
(largely due to the Majorana phase factor /) of the allowed bands in the two lower panels. 
It will then be clear that probing the spectrum hierarchy and the Majorana phase(s) with 
future experiments will be a very challenging task. 



VIII. ADDING LOWER BOUNDS ON m^p FROM THE CLAIMED Oi/2/3 SIGNAL 

In this Section we keep the same data set as in Sec. VII, except for the Oi/2/3 data input, 
which we now take from Eq. in order to show the effect of the signal claimed in 4^ 3 
on the global fit. 

Figure 5 shows that, in this case, there is a lower bound on mpj^ at 2cr, as indicated by 
an additional horizontal dashed line (not present in the previous Fig. 4). This lower bound 
is somewhat "too high" to allow a good combined fit with oscillation and cosmological 
CMB+LSS data (beta decay limits still being not relevant in the global fit). Therefore, the 
global 2(7 allowed region extends somewhat outside the 2a limits from cosmology and Oj/2/3 
data taken separately — a clear indication of some statistical tension between such data. 

From Fig. 5 it follows that, assuming standard 3z/ mixing and the 0i/2/? claim in 4^ 4^ 1 : 



1) The preferred spectrum is degenerate and there is no possibility to discriminate normal 
and inverted hierarchy in the (771/3,771/3/3, S) parameter space; 2) A significant fraction of the 
2(7 region allowed by the current global fit is within the sensitivity reach (771/3 ~ 0.2-0.3 eV) 
of the future beta-decay experiment KATRIN |38>]; 3) A cosmological signal for S ~ 0(1) eV 
should be "around the corner;" 4) if such a signal will not be found, the tension between 
Oi/2/3 and S data will increase, unless the theoretical central value for the matrix element 
Cmm is significantly shifted upwards, so as to decrease the preferred values of 711/3/3. Indeed, 
we shall see in the next section that such tension is already increased by including recent 
Lya data. 

Let us speculate, however, about the possibility that the true values of {mp^mpf^^Ti) lye 
within the allowed regions of Fig. 5. It appears that, in order to significantly reduce such 
regions, one should: improve the sensitivity to mp by a factor ~ 5 at least, find a cosmological 
signal of S ~ 0(eV) with an uncertainty definitely smaller than a factor ~ 2, and confirm the 
current Oz/2/5 signal claim with a total (experimental+theoretical) error reduced by better 
than a factor of 2 (with respect to our estimate). Improvements in neutrino oscillation data 
would instead produce marginal effects (not shown). These stringent requirements set the 
stage for possible future tests of the absolute neutrino mass scenario allowed by the global 
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data fit in Fig. 5. 



IX. ADDING Lya FOREST DATA 

In this Section we re-evaluate the bounds shown in Sec. VII (without 0y2d signal) and 
in Sec. VIII (with Oz/2/3 signal), by adding the recent Lya forest data jsj IH^ discussed in 
Sec. V. 

A. Case with no 0u2(3 signal 

Figure 6 shows the global fit assuming no 0^2(3 signal, but adding Lya data. The same 
fit, without Lya data, has previously shown in Fig. 4. By comparing Figs. 4 and 6 the 
impact of Lya data, taken at face value, is impressive: The upper limit on E is improved by 
a factor ~ 3 and, through the correlations induced by neutrino oscillation data constraints, it 
is transformed into upper limits onto and m^^, which are an order of magnitude stronger 
than the current laboratory ones. The overall bounds are strong enough to approach the 
regime of partially degenerate spectrum, where the bands for normal and inverted hierarchies 
start to branch out in the two left panels of Fig. 6. 

The bounds in Fig. 6 set very stringent requirements to future laboratory experiments 
sensitive to absolute neutrino mass. A factor of ^ 10 improvement is required both for rrif^ 
and for measurements, in order to be competitive with the cosmological upper bounds 
and to improve current upper limits in the {mp^mpf^.T^) parameter space; finding a signal 
would be even more demanding. On the other hand, cosmological data should eventually 
provide evidence for nonzero S in the range 0.05 ^ S/eV ^ 0.5 for normal hierarchy, or in 
the range 0.09 ^ S/eV ^ 0.5 for inverted hierarchy. 

B. Case with claimed Oz^2/3 signal 

The strong upper bound on S obtained by adding Lya data to the cosmological data 
input increases the already existing tension with the Oz/2/5 claim (see Sec. VIII) to the point 
that a global fit would provide meaningless and unstable results. 

For this reason, we show in Figure 7 the bands separately allowed at 2cr by the 
Oz/2/3 claim (horizontal band) and by the combination of all oscillation and cosmological 
(CMB+2dF-|-Lya) neutrino data (lower slanted bands), for both normal and inverted mass 
hierarchy. Only the plane {mpp, E) is shown in Fig. 7, since current laboratory bounds on 
171/3 are an order of magnitude away from the global allowed region, as previously discussed. 

In Fig. 7, the absence of overlap between the bands separately allowed at 2a is a clear 
symptom of possible problems, either in some data sets or in their theoretical interpretation, 
which definitely prevent any global combination of data. However, it would be premature to 
conclude that, e.g., the 0^2^ claim is "ruled out" by cosmological data. Firstly, cosmological 
bounds on S are rather indirect, and involve the processing of a large amount of data, 
whose systematics must be more closely scrutinized, in particular for the latest Lya data 
set. Secondly, one cannot exclude that future calculations of the nuclear matrix element 
Cmm in Oz/2/5 may be revised upwards, thus lowering the m/3/3 allowed range and relaxing 
the tension with current upper limits on S. Thirdly, one cannot exclude that Oi^2/5 decay 
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might receive contributions from new physics effects beyond the exchange of hght Majorana 
neutrinos. Finally, either some data or some fundamental assumptions about the standard 
three-neutrino mixing and the cosmological scenarios might be wrong. With all the above 
cautionary remarks, it is anyway exciting that global neutrino data analyses have already 
reached a point where such fundamental questions start to arise. 

X. CONCLUSIONS AND PERSPECTIVES 

In the context of standard 3z/ mixing and standard cosmology, we have performed a 
global phenomenological analysis of the constraints applicable to three observables sensitive 
to absolute neutrino masses: the effective neutrino mass in Tritium beta decay (mp); the 
effective Majorana neutrino mass in neutrinoless double beta decay {mpp)] and the sum of 
neutrino masses in cosmology (S). 

We have first discussed the correlations among such variables induced by neutrino os- 
cillation data, in both normal and inverse neutrino mass hierarchy (see Fig. 3 and related 
comments). We have then applied laboratory constraints on mp and m^j/j, as well as cosmo- 
logical constraints on E, in order to obtain global fits in the (m^, E) parameter space, 
which embed all world neutrino data relevant to absohitc neutrino mass scenarios. The 
results have been discussed in terms of two-dimensional projections of the globally allowed 
region in the {mp^mpp, E) parameter space, which neatly show the relative impact of each 
data set. In particular, the (in) compatibility between E and mgQ constraints has been dis- 
cussed for various combinations of data, as shown in Figs. 4 and 5 (without and with the 
Qv2j3 signal claim, respectively) and in Figs. 6 and 7 (which include recent Lya data). 

We have also briefly discussed how future neutrino data (both oscillatory and non- 
oscillatory) with improved sensitivity can further probe the currently allowed regions. Our 
graphical representations appear to be rather useful in this sense, since prospective measure- 
ments of any of the three observable (m^, m^^, E) can be easily mapped onto the currently 
allowed regions. 
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3]y oscillation parameter bounds fronn marglnolized A^o: 




FIG. 1: Bounds on the Su mass-mixing neutrino parameters ((5m^, Am^, sin^ ^12, sin^ ^13) from our 
analysis of current neutrino oscillation data. The solid curves represent projections of the Axosc 
fit function. The intercepts of the curves with the horizontal dotted lines at Ax^ = v? define the 
n-a limits on each parameter. There is no statistically significant lower limit for sin^ ^13. 
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FIG. 2: Upper bounds on the sum of neutrino masses S from our 3z/ analysis of cosmological 
data, given in terms of the A^l function. The sohd and dashed curves refer to the combination of 
CMB and LSS data (CMB+2dF and CMB+SDSS, respectively). The two CMB+LSS fits provide 
comparable results and, for definiteness, the CMB+2df one is adopted. In addition, we consider also 
the case where the recent Lya data from the SDSS are included, providing significantly stronger 
constraints on S (dotted curve). See the text for details. 
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2a bounds from I 
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FIG. 3: Global 3z/ analysis in the {mp^mj^p^Ti) parameter space, using oscillation data only. The 
three panels show the two-dimensional projections of the volume allowed at A^osc = 4 (2a on 
each parameter) for both normal hierarchy (thick solid curves) and inverted hierarchy (thin solid 
curves), respectively. The emerging correlations between parameters are discussed in the text. 
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FIG. 4: Global 3z/ analysis in the {mp^nipp^Ti) parameter space, using oscillation data plus 
laboratory data and cosmological data. With respect to Fig. 3, this figure implements also upper 
limits (shown as dashed lines at 2a level) on from Mainz+Troitsk data, on mpp from Oi^2/3 
data, and on S from CMB+2dF data. In combination with oscillation parameter bounds, the 
cosmological upper limit on E dominates over the laboratory upper limits on mp and m^p. 
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FIG. 5: Global 3z/ analysis in the (m^,m^^,S) parameter space, using oscillation data plus 
laboratory and cosmological data. With respect to Fig. 4, lower bounds on m^p from a claimed 
Qv2(5 signal are implemented. The globally allowed 2a regions appear to be stretched somewhat 
beyond the separate 2a bands. This feature reflects some tension existing between the Oi/2/3 claim 
and cosmological CMB+LSS data within the oscillation scenario. 
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FIG. 6: Global 2>v analysis in the {mp,mpp,Ti) parameter space as in Fig. 4, but including Lya 
forest data from the SDSS, which lead to a stronger upper bound on S, close to the zone where 
normal and inverted hierarchies start to branch out in the (m^,S) and (m^^,S) planes. 
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FIG. 7: 3u analysis in the (m^^, S) plane, showing the absence of overlap between the regions 
separately allowed at 2a by the Oi^2/3 claim (horizontal band) and by the combination of u oscillation 
data with CMB+2dF+Lya cosmological data (slanted bands), for both normal hierarchy (N.H.) 
and inverted hierarchy (I.H.). Such strong tension might indicate possible problems either in some 
experimental data or in their theoretical interpretation. 
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